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Abstract — A microelectronic circuit of block-elements 
functionally analogous to two hydrogen bonding networks is 
investigated. The hydrogen bonding networks are extracted 
from a-lactamase protein and are formed in its active site. 
Each hydrogen bond of the network is described in equivalent 
electrical circuit by three or four-terminal block-element. 
Each block-element is coded in Matlab. Static and dynamic 
analyses are performed. The resultant microelectronic circuit 
analogous to the hydrogen bonding network operates as 
current mirror, sine pulse source, triangular pulse source as 
well as signal modulator. 

Index Terms — Hydrogen bonding network, behavioral 
modeling, Matlab, Verilog-A, proteins. 

I. Introduction 

Bioelectronics is dynamically evolving part of electronics. 
The integration of biomolecules with electronic elements to 
yield functional devices attracts substantial research efforts 
because of the basic fundamental scientific questions and 
the potential practical applications of the systems. A 
fundamental requirement of any bioelectronics system is the 
existence of electronic coupling and communication between 
the biomolecules and the electronic supports [1]. The 
understanding of charge transport phenomena through 
biological matrices attracted in the past decades, and 
continues to evolve, intensive theoretical and experimental 
work. The seminal contributions of the Marcus theory [2], 
the super exchange charge transfer theory [3], and the 
definition of superior tunneling paths in proteins [4] had a 
tremendous impact on the understanding of biological 
processes such as the electron transfer in the photosynthetic 
reaction center, or the charge transport in redox -proteins that 
are the key reactions for numerous electrochemical and 
photoelectrochemical biosensing systems. 

The small dimensions of the bioobjects make them 
appropriate for various types of biosensors [5], pH-sensors 
[6], light sensors [7], etc. A biosensor is generally defined as 
an analytical device which converts a biological response 
into a quantifiable and processable signal. A continuous 
feedback between elegant experimental work employing 
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structurally engineered proteins and theoretical analysis of 
the results led to the formulation of a comprehensive paradigm 
for electron and proton transport in proteins [8]. 

In our previous research [9, 10, 11, 12, 13] we have 
investigated information transfer through hydrogen bonding 
networks of the /^lactamase protein. After extraction of protein 
hydrogen bonding networks (HBNs) from the enzyme 
periphery, we have investigated the charge transport. We 
have matched each residue (from protein HBN) to functionally 
analogous block-element to describe the proton transfer with 
polynomials. Next, static and dynamic analyses have been 
performed on the circuits consisting of block-elements. The 
simulation results are compared to conventional 
microelectronic circuits. 

In this paper we focus on information transfer via protein 
hydrogen bonding networks and their operation as analogs 
to conventional electronic circuits. In particular, we 
investigate HBNs formed in the active site of the ^-lactamase 
protein. 

II. Model and Equations 

The two hydrogen bonding networks formed at the ac- 
tive site during the four intermediates of acyl-enzyme reac- 
tion are shown in Fig. 1 ; they are denoted "nucleophilic" and 
"electrophilic". Proton transfer in each hydrogen bonding 
network was studied in [14]. In the paper is introduced pro- 
ton transfer parameter K (in accordance with the constant of 
the electron transfer). It has dimension of free energy, and 
from the calculations it can be interpreted as follows: the 
greater parameter K - so much readily accomplished proton 
transfer between donor and acceptor from HBNs, i.e. the pro- 
ton current will be greater. On the other hand the parameter 
of proton transfer depends on the donor/acceptor potentials 
similarly to the potentials supplying the microelectronic com- 
ponents. Therefore we can construct three and four-terminal 
electronic block-elements analogous to the hydrogen bonds. 

During the four intermediates of the half-cycle of the acyl- 
enzyme reaction participate the same residues and water 
molecules. For this reason the corresponding microelectronic 
block-elements in the equivalent electrical circuits are the 
same but they are connected in different way. Therefore, two 
functionally analogous microelectronic circuits for the four 
intermediates of acyl-reaction are constructed; the circuits 
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Figure 1. Acyl-enzyme reaction intermediates - HBNs in the active site of: (E) free enzyme, (ES) Michaelis complex, (Tl) transient state, 
(EY) acylenzyme. The first HBN is referred to as "nucleophilic" consists of residue S70, water molecule w297, residues N170, E166, K173, 
N132 and ligand. The second referred to as "electrophilic" consists of consists of residues S130, K234, water molecule w309, and residues 

D214, S235 



depend on switch S that corresponds to the ligand. In this 
paper we will analyze one intermediate of the acyl-reaction 
only, i.e. the free enzyme (state (E) in Fig.l) without a ligand. 

In Fig.2 it is shown the microelectronic circuit analogous 
to HBNs formed in the active site of p-lactamase. 
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Figure 2. Microelectronic circuit analogous to HBNs formed in the 
active site of /^-lactamase when /^-lactamase is free enzyme 
(state E) 

The input voltage of the first electric circuit analogous to 
"nucleophilic" hydrogen bonding network is labeled as 
"Uinl". Block-element Tl corresponds to lysine K73NZ. The 
lysine is a strong proton donor, so in the circuit it is interpreted 
as a voltage controlled current source. Block-element T2 is 
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juxtaposed to S70OG residue, which is a three terminal block- 
element with identical input and output voltage and different 
input and output currents. 

In the circuit of Fig. 2 the water molecule w297 is compared 
to block-element T3. It should be noted that depending on 
the state of the switch S, block-element T3 changes its 
function. In this case T3 has two inputs and one output. 

The block-element T4 is functionally analogous to El 66, 
which is strong proton acceptor and can form different 
hydrogen bonds. T4 block-element is represented as a three 
terminal block-element with one input and one output. The 
other two block-elements T5 and T6 are similar to asparagines 
N132 and N170. Asparagine could be both proton-donor and 
proton-acceptors in the circuit, therefore N132 and N170 can 
change the direction of the currents as a depending on the 
switch S. 

The input of the second circuit on Fig. 2 is labeled with 
Uin2. This electrical circuit is functionally analogous to the 
"electrophilic" hydrogen bonding network, which also has 
three outputs. Here block-element T10 is also a voltage 
controlled current source and it corresponds to the strong 
proton-donor K234NZ. In this block-element (similarly to Tl) 
the two output voltages are identical to the input voltage but 
the input and output currents are different. The next block- 
element Til is analogous to the water molecule w309. In 
state (E) the element has two inputs and one output. The 
residue D214 is juxtaposed to block-elements T12; T12 has 
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one input and one output, and the T12 output is also the 
circuit output. The three terminal block-element T13 is 
compared to S130OG. Its input and output voltages are 
identical, but the input-output currents are different. The 
block-element T14 is the last output block-element of the 
integrated circuit. It is functionally analogous to S2350G and 
has the same functions as S130OG (respectively T 13). 

It should be noted that there is proton transport in both 
HBNs during the different intermediates of the reaction. The 
proton transport is influenced by pH and protein-ligand 
interaction. In this reason, both electrical circuits are not 
separated from each other (similar to an integrated circuit) 
and the input voltages of both units are equal in analogy of 
pH. 

The model equations that describe the electrical relations 
of the each block-element are given below. 

Equations (1) and (2) describe the Tl block-element. 



77, = 7/ 

1 in 



(1) 



7, = 3x10 5 7// - 5x10 - 5 Uf -7x10 5 77, + 0.0013; (2) 



Equations (3) and (4) describe the T2 block-element. 



U 2 = 1.0451x77,-0.1194; 



(3) 



/„ = 2xl0 5 x7/ n 3 -13xl0 6 x7/ n 2 -0.00015x7/,, 

+0.00048; (16) 
Equations (17) and (18) describe the T12 block-element. 

7/ =0.9683x7/. , +0.458; (17) 



(18) 



Equations (19) and (20) describe the T13 block-element. 

7/, 3 = 1.1009x7/, -0.3571; (19) 



/ 13 = 10 5 x7/, 3 4 +2xl0 5 x7/, 3 3 -8xl0 5 7/,/+2xl0- 5 x 7/, 3 



+0.0015; (20) 
Equations (21) and (22) describe the T14 block-element. 

U= 1.034xt/„ -0.1986; (21) 



7 14 =0x77 M +0.0001; 



(22) 



The equations are coded in Matlab [15]. Excerpt of the 
code is given below. 



7, =-0.0236x7/ 3 + 0.0379x7// + 0.124x7/, +1.3046; (4) 



Equations (5) and (6) describe the T3 block-element. 

7/ 3 = 1.0179x7/ 2 - 0.039; (5) 

7 3 = 0.0015x77/ - 0.0015x77/ -0.0171xT/ 3 + 0.0859; (6) 
Equations (7) and (8) describe block-element T4. 

77 4 = 0.9703x7/ 3 +0.0589; (7) 



4 3 outl' 



(8) 



Equations (9) and (10) describe the T5 block-element. 

U 5 = 0.0457x77/ +1.2273x77,-0.8501; (9) 



7 5 = Ox U s +0.0001; 



(10) 



Equations (11) and (12) describe the T6 block-element. 

(11) 



77 =1.0544x77 -0.0933; 



7 6 = 0.0058x77/ + 0.0367x7/ 6 + 0.71 13; (12) 
Equations (13) and (14) describe the T10 block-element. 



77, =0.9701x77, +0.2072; 



7 ]0 = 4x10 5 x7/, 4 - 6x10 5 x7/,/ - 0.0002x7/,/ + 



0.0004xT7, n + 0.00252; 



(13) 



(14) 



Equations (15) and (16) describe the Tl 1 block-element. 

(15) 



7/ n = 0.9835x7/ I0 +0.1438; 



U2=1.0451*U1 -0.1194; 

plot(Ul ,U2, 'linewidth' ,2); 

set(gca, ' fontweight' , ' b' , ' fontsize' ,14) 

grid on 

title('T2'); 

xlabel('Uin[V]'); 

ylabel('U2[V]'); 

legend( ' simulation' , ' data' ) ; 

set(legend( ' simulation' , ' data' , 1 ), ' fontsize' , 12) ; 

pause; 

12 = - 0.0236*U2. A 3 + 0.0379*U2. A 2 + 0. 124*U2 + 1 .3046; 
load('bsyl_e.dat'); 
U2exp = bsyl_e(:,3); 
I2exp = bsyl_e(:,4); 

plot(U2,I2,U2exp,I2exp,'ro','linewidth',2); 

set(gca, ' fontweight' , ' b' , ' fontsize' , 14) 

grid on 

title('T2'); 

xlabel('U2[V]'); 

ylabel(T2[pA]'); 

% legend('simulation','data'); 

set(legend( ' simulation' , ' data' , 1 ), ' fontsize' ,12); 

pause; 

III. Static Analysis 

Static analysis is performed in Matlab. First, we tested 
how well our polynomial 7- V characteristics describe the com- 
ponents and connections between them. In Fig. 3 are shown 
the compared results from simulations and the previous data 
for block elements Tl and T2. 

As it can be seen from Fig. 3, the polynomials well describe 
the functional dependencies of the modeled block elements. 
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Figure 3. I-V characteristic of block elements a) 01 and b) 02 
(emulating El 3, S70) 

Similar results were seen for the other block elements of the 
circuit (by this reason they are not shown), and the maximum 
calculated error is 5.66%. 

After the comparative analysis between the simulations 
and the data taken from [14], a static analysis is made for the 
entire integrated biocircuit. The results for each output of 
the circuit are shown below. Fig. 4 shows the /- V characteris- 
tics of output 1 . 
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Figure 4. I-V characteristic of first output current 
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The I-V characteristic of first output is similar to I-V 
characteristic of tunnel diode [9]. 

In Fig. 5 the output voltages are shown. 



Uouts vs. Uin 




Figure 5. The output voltages vs. input voltage. 

The output voltages repeat the form of the input voltage. 
This is also observed in our previous paper for HBNs in the 
periphery of /^-lactamase protein [9] . 

In Fig. 6 the current of the second output is shown. 
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Figure 6. I-V characteristic of second output current. 

Current does not change with the change of voltage (its 
value is 0.0001 [pA]). Therefore, this output has 
characteristics similar to a current source as observed in [10, 
11]. 

In Fig. 7 the current of the third output is shown. 

For case E current increases according to an exponential 
law with increasing voltage and output characteristic is similar 
to shifted characteristics of diode. 

Fig. 8 shows the current characteristic on output 12. 

Current is amended by S-shaped law. 

At the next output (13) there is a curve which is a 
combination of the curves of S-shape (Fig. 9). 

At the next output 14 (Fig. 10) the current does not 
change by changing the voltage. 

Characteristics are similar to the current generator (Fig. 
10), as it was for output 2 [12]. 

From the standpoint of its mode of use, the circuit could 
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Figure 9. I-V characteristic of thirteenth output current 

be used as a current generator or as an amplifier. Further its 
capabilities will be demonstrated in the dynamic analysis. 

IV. Dynamic Analysis 

The integrated circuit functionally similar to hydrogen 
bonding networks in the active site of the protein is a complex 
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Figure 10. I-V characteristic of fourteenth output current 

circuit with multiple output; therefore, a thorough dynamic 
analysis is required. To understand in detail its mode of action, 
research were made when applying three types of sinusoidal 
voltages at the input. First was applied voltage with an 
amplitude from -2.2 to 2.2 [V] (Fig. 11) and the frequency of 10 
[GHz] then was applied voltage with the same frequency, but 
only the positive or only the negative amplitudes and sweep 
of 1 [V] . These simulations showed that the output voltage 
repeats the input voltage in form and frequency, and the 
amplitude is on the same order. 
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Figurell. The input voltage vs. time Uin = 2.2*sin((5el 1 )*t) 

This behavior was expected because from static analysis 
we know that the output voltages are shifted off the input 
voltage by linear law. Therefore, for convenience we will 
present only the input voltage and output currents vs. time. 
For output currents again we see that they keep their posi- 
tive amplitude regardless whether voltages are in negative or 
positive half-period. As an example, in Fig. 12 we showed 
voltage and current on output 1, input voltage is with ampli- 
tude from -2.2 to +2.2 [V]. If input voltage with positive ampli- 
tude is applied (Fig. 13, 14) different signal amplitude sinu- 
soidal signals are generated at output 1; the currents are 
shifted in relation to the voltages. 



97 



—ACEEE 



Short Paper 



ACEEEInt. J. on Signal and Image Processing , Vol. 5, No. 1, January 2014 



T4 OUT1 



T4 OUT1 



3 
2 
1 

s 

5 

3 
O 

3 

-1 











-simulation 




















































I 0.2 


0.4 


0.6 


0.8 1 



t[s] 
T4 OUT1 



X10 



0.11 
0.1 
0.09 
0.08 
0.07 
0.06 
0.05 













simulation 











































































0.2 



0.4 0.6 
t[s] 



0.8 



x10 



Figurel2. The voltage and current on first output vs. time (Uin : 
2.2*sin((5el l)*t 
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Figure 13. The input voltage vs. time Uin =1,5 +0.5 *sin((5el l)*t) 

When applied a input voltage with negative amplitude 
(Fig. 15) for output 1, we obtain the modulated signal with 
the same phase for current (Fig. 16). Similar phenomena was 
observed in our previous paper [13]. 

Simulations show that the second output current does 
not change (Fig. 17) when applying with different amplitude 
(and frequency) input voltages. 

This output is similar to the output of the current-mirror, 
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Figure 14. The first output current vs. time (Uin =1,5 
+0.5*sin((5el l)*t)) 
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Figure 15. The input voltage vs. time Uin = -l+0.5sin((5el l)*t) 
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Figure 16. The first output current vs. time Uin = - 
l+0.5sin((5el l)*t) 

i.e. the currents not dependent on a change of the voltage, 
as determined by the static analysis. 

When applying a sinusoidal voltage with a positive am- 
plitude (Fig. 15) the output 3 starts to generate sinusoidal 
signals (Fig. 18). 

The same output signals are obtained when applying the 
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Figure 17. The second output current vs. time 
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Figure 18. The third output current vs. time (Uin 
+0.5sin((5ell)*t)) 
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= 1,5 



sinusoidal input voltage with negative amplitude. Therefore 
they are not shown in any figures. 

The most complicated generated signals are obtained at 
the output 12. When is applied a sinusoidal input signal with 
an amplitude of -2.2 to +2,2 V, then output starts to modulate 
(Fig. 19). 

When applying an input voltages with only positive or 
only negative amplitude, the circuit allows the generation of 
other types of signals (triangular) as shown below (Fig. 20). 

Generally, only this output of the circuit is able to generate 
signals, such as number of previously known in electronic 
devices. These signals can be controlled by changing the 
amplitude of the input signal. 

The situation is similar at output 1 3 and there is generation 
of signals with various shape and amplitude depending on 
the amplitude of the input voltage (Fig. 21, 22). The 
comparison of the signals at output 13 and output 12 showed 
that these signals are different. 

Output T13 can easily be put into mode that generates 
sinusoidal signals only. This happens when applying the 
input voltage with a negative amplitude (Fig. 17). 

Here in Fig. 22 current is in phase with the input voltage 
(respectively the output voltages). 

The output current at T14 output does not change with 
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Figure 19. The output current on output 12 vs. time (Uin 
2.2*sin((5ell)*t)) 
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Figure 20. Other types of signals generated on output 12. 

the input voltage (Fig. 23): the currents are constant (0.0001 
[pA]) and the output behaves similarly to the current-mirror 
or current source. 

The results of the dynamic analysis showed that the 
integrated circuit functionally analogous to the hydrogen 
bonding network in the active site of the protein /^lactamase 
could operate in different modes and to generate different 
types of signals and can also modulate the different in 
amplitude and shape signals. 
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Figure 21. The output current on output 13 vs. time (Uin 
2.2*sin((5ell)*t)) 
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Figure 22. The input voltage and output current on output 13 vs. 
time (Uin = -1.5 +0.5sin((5el 1 )*t)) 
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Figure 23. The output current on output 14 vs. time (Uin = 
2.2*sin((5ell)*t)) 

V. Conclusion 

The block-elements of the microelectronic circuit well 
describe the behavior of the hydrogen bonding network. Static 
analysis shows that the circuit could operate as tunnel diode, 
current source and multifunctional device. From the results 
of the dynamic analysis, it can be seen that the circuit could 
operate in different modes and generate in their outputs 
different types of signals. The circuit has functions similar to 
current-mirror, source of different signals with sinusoidal and 
triangular amplitude, and can modulate signals with different 
amplitude and shape. The results obtained prove that the 
modeled circuit can emulate the behavior of hydrogen 
bonding networks for microelectronics applications. 
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